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MINOS will measure the composition of a beam at two locations, 735 km apart, in an effort confirm the 
(atmospheric) neutrino oscillation hypothesis and measure the mixing parameters Am^ and sin 2 (26*23). Oscilla- 
tions will be manifested as a difference in the rate and energy spectrum of — CC interactions measured in the 
two detectors. Because most interactions observed in MINOS are inelastic, the f M energy is reconstructed as the 
sum of the energy carried by the muon and that seen in the hadronic shower emanating from the struck nucleus. 
The latter is sensitive to uncertainties in the hadronisation process, chief among them those due to intranuclear 
re-scattering (i.e., final state interactions). We discuss the simulation of intranuclear re-scattering currently used 
by MINOS and its effect on quantities observable in the experiment. 



1. Introduction 

The Main Injector Neutrino Oscillation Search 
(MINOS) is a long baseline, two-detector neu- 
trino oscillation experiment that will use a muon 
neutrino beam produced by the Neutrinos at the 
Main Injector (NuMI) facility at Fermi National 
Accelerator Laboratory (FNAL) P21- The mea- 
surement will be conducted by two functionally 
identical detectors, located at two sites, the Near 
Detector (ND) at FNAL and the Far Detector 
(FD) in the Soudan Underground Laboratory in 
Minnesota. The NuMI beamline j5] and the 
735 km long-baseline will allow exploration of the 
A7TI23 and sin 2 (2^23) mixing parameters studied 
previously with atmospheric neutrinos 0] and by 
the K2K experiment [5]. For — > v e transitions, 
MINOS will probe neutrino mixing parameters 
beyond the current limits of the CHOOZ experi- 
ment 0. 

The MINOS detectors are tracking-sampling 
calorimeters, optimised to measure neutrino in- 
teractions in the energy range 1 < E v < 
50GeV. The active medium comprises 4.1 cm- 
wide, 1.0 cm-thick plastic scintillator strips ar- 
ranged side by side into planes. Each scintillator 
plane is encased within aluminum sheets to form 
a light-tight module and then mounted on a steel 
absorber plate. The detectors are composed of a 



series of these steel-scintillator planes hung ver- 
tically at a 5.94cm pitch with successive planes 
rotated by 90° to measure the three dimensional 
event topology. Wavelength-shifting and clear op- 
tical fibers transport scintillation light from each 
strip to Hamamatsu multi-anode photomultiplier 
tubes which reside in light-tight boxes alongside 
the detector. Both detectors are magnetized so 
as to measure muon charge-sign and momentum 
via curvature. 

The disappearance measurement will be 
done by using the event topology to identify in- 
teractions as — CC, rather than NC/^ C — CC, 
then reconstructing the neutrino energy as the 
sum of the muon energy and shower energy: E u = 
Eft + Shad- The oscillation hypothesis will be 
tested by comparing the measured neutrino en- 
ergy spectrum to the spectrum expected in the 
absence of oscillations. The latter spectrum will 
be anchored to observations made with the Near 
detector. MINOS endeavors, in the case of oscil- 
lations, to measure Am| 3 and sin 2 (2^23) with an 
accuracy of better than 10%. The results depend 
upon a reliable knowledge of the event selection 
efficiency and the energy scale, both of which may 
be influenced by uncertainties in the hadronisa- 
tion process. 
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Figure 1. Two simulated events in the MINOS Near Detector, zoomed to show the vertex region. The 
large arrow marks the direction of the incoming neutrino, lines mark the direction of the outgoing muon 
(blue), protons (red) and pions (pink), and dots show hits recorded by the detector. The dark vertical 
bands denote the iron absorber plates (not to scale). MINOS measures two orthogonal "views" of each 
event of which only one is shown here. 



2. Monte Carlo Simulation 

We will show results from a full Monte Carlo 
simulation of the MINOS experiment. The sim- 
ulation consists of four parts: (1) a neutrino tar- 
get and beamline simulation (GNuMl), (2) the neu- 
trino interaction generator NEUGEN [7], (3) a full 
detector GEANT-based physics simulation known 
as GMINOS and (4) a detailed active detector re- 
sponse simulation. The detector physics and re- 
sponse simulation has been benchmarked against 
cosmic ray and beam muons measured in the Near 
and Far detectors as well as series of testbeam 
measurements collected with a scaled down ver- 
sion of the MINOS detectors in the CERN PS 
East Hall. The test beam measurements are used 
to fix the energy scale [Hj (via muon range) and 
validate the simulation of electromagnetic JSJ and 
hadronic showers |10| . The hadronic shower code 
GCALOR was found to be in good agreement with 
the data and is used in our simulations. 

Neutrino interactions are modelled by 



NEUGEN- v3. NEUGEN was originally written to 
simulate neutrino backgrounds in the Soudan-2 
proton decay experiment. Since that time it has 
undergone substantial modification to upgrade it 
for use at higher energies and has been adopted 
by MINOS. NEUGEN includes a Rein-Seghal based 
treatment of neutrino induced resonance produc- 
tion, CC and NC coherent pion production and 
a modified leading order DIS model extended 
to improve the treatment in the transition re- 
gion between DIS and resonant production 
KNO scaling is used to model the final state 
multiplicity in the DIS regime. 

Two simulated events in the are shown in Fig^ 
The detector has a relatively fine segmentation for 
a neutrino calorimeter but was optimised to max- 
imise the event rate in the FD, equalise perfor- 
mance of the two detectors, discriminate between 
v/v, and cover a wide range of neutrino energies. 
The detector is not able to measure details of the 
final state in most interactions. Figure Hal is no- 
table in that it contains a short but distinctive 
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track-stub corresponding to a ~ 1 GeV/c proton. 
This proton is near the lower practical limit for 
identification in quasi-elastic events. 

NEUGEN also includes an intranuclear re- 
scattering package known as INTRANUKE. The 
code is anchored to a comparison of final states in 
v + d and v + Ne interactions as measured in the 
BEBC and ANL-12ft bubble chambers P2J. The 
library includes a treatment of pion elastic and 
inelastic scattering, single charge exchange and 
absorption in a cascade simulation of the final 
state. The relative probabilities for these pro- 
cesses are approximately 35:50:10:5 (40:30:7:23) 
at a pion kinetic energy of 1 GeV (250 MeV). 
The formation zone concept is included by sup- 
pressing interactions that would have occurred 
before the pion has travelled a distance I = 
Top/m To = 0.52 fm/c |13I14) . The simulation 
includes a treatment of pion absorption inspired 
by |15j . Absorbed pions transfer their energy to a 
npnp cluster and the individual nucleons are then 
tracked in GEANT. 

3. Results 

Figure [21 shows the probability for absorption 
and inelastic scattering as a function of the ini- 
tial pion momentum. Some familiar features, 
such as the A peak can be clearly seen. The 
decrease in the scattering probability with pion 
energy is a consequence of the formation zone. 
The effective absorption probability is signifi- 
cantly larger than the relative probability stated 
above. This is caused by the fact that the as- 
sumed nuclear radius is much larger than than 
the mean free path, making multiple interactions 
common. Also, more importantly, pion absorp- 
tion is a "one way street" in the sense that the 
cascade ends and other processes are suppressed. 

Pions lose a large fraction of their momentum 
in the intranuclear re-scattering process. Figure^ 
shows the average loss as a function of the pion 
momentum for all pions exiting the struck nu- 
cleus and separately for only those pions that un- 
derwent an inelastic interaction. The underly- 
ing distribution in the latter case is quite broad. 
The simulation assumes that energy lost in in- 
elastic interactions is predominantly transferred 



to heavy remnants of the target nucleus which 
are expected to range out in the iron absorber, 
remaining invisible to the experiment. 

It is clear that low energy, fully active, neu- 
trino scattering experiments seeking to measure 
exclusive final states will be rather sensitive to re- 
scattering. However, based on the rather high in- 
teraction probability and significant energy loss, 
one should also expect observable consequences 
even in a dense neutrino calorimeter like MINOS. 
This is indeed the case. Figure0]shows the change 
in the measured vertex shower energy as a func- 
tion of the energy transfer v. The origin of the 
vertical axis corresponds to the default configura- 
tion employed in our simulations. We observe a 
substantial change in the measured shower energy 
under different re-scattering scenarios. 

The difference between the default simulation 
and that without any re-scattering indicates that 
as much as 12% of the initial shower is expected to 
go unobserved. Neutrino oscillation experiments 
that directly use E v depend upon reconstruct- 
ing the energy (see Fig 03) in an unbiased fashion 
or, at the very least, applying a correction fac- 
tor when extracting the oscillation parameters. 
In this regard MINOS is in a somewhat better 
position than its predecessors since the measured 
spectrum in the Near Detector provides a strong 
constraint when testing for oscillations and ex- 
tracting the mixing parameters. However utilis- 
ing this data requires a reasonable understanding 
of the accuracy of the rescattering formulation in 
order to disentagle its effects from those caused 
by uncertainties in the neutrino flux and cross- 
sections. 

In Fig 2| "Full 7r absorption" corresponds to ab- 
sorption with the same probability as in Fig [2] 
but in which the 7r energy is not redistributed to 
an npnp cluster. It is clear that the details of 
the pion absorption process are very important 
and that the shower energy measurement is sen- 
sitive to nucleons with a kinetic energy of only a 
few hundred MeV. Figure 0] also shows that the 
effect of final state interactions is not confined 
only to low shower energies. This is because even 
in large hadronic showers the energy is carried 
by multiple, lower energy particles. We have ex- 
plored this effect by simulating with formation 
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times r = 0.34 fm/c and 0.80 fm/c. This is a 
large variation but only results in a 2-3% modifi- 
cation of the detector response. 

Since MINOS must use the event topology to 
remove the NC background we have studied the 
effect that re-scattering has on the efficiency and 
purity of the sample. The selection procedure 
chooses events as v^ — CC based on the presence 
of a track with a length consistent with a muon 
having a momentum larger than 500 MeV/c, one 
or or zero showers near the track vertex, a good fit 
to the curvature in the magnetic field and a track 
pulse-height consistent with one MIP per plane 
crossing. We find that with this selection the 
i^— CC efficiency and the level of NC background 
varies no more than a 1-2% over the energy range 
1 < E v < 20 GcV. Such small variations can be 
constrained using data from the Near detector. 

Finally, a more interesting situation occurs 
when we attempt to select a quasi-enriched sub- 
sample of events according to a procedure that 
uses the reconstructed invariant mass, the pulse- 
height not assigned to the track, the fraction of 
hits with a pulse-height larger than 20 PE and 
the number of hits assigned to secondary tracks 
found via a Hough transform. Because baryon 
interactions are not simulated, the efficiency of 
this selection is independent of the final state in- 
teraction scenario. The purity of the selection, 
however, is not. Without re-scattering the purity 
of our sample integrated over 0-10 GeV was 80%, 
under the n — > npnp scenario it decreased to 73% 
and with full n absorption only 63% . This be- 
haviour is broadly in line with expectations as our 
procedure tends to select events with low vertex 
activity as quasi-elastic. 

4. Conclusions 

We have explored the ramifications of intranu- 
clear re-scattering in MINOS through Monte 
Carlo studies using the NEUGEN event generator 
coupled to a GEANT based detector simulation. 
The results of the study indicate that the re- 
constructed shower energy and topology depends 
rather strongly on the amount of energy lost in 
the re-scattering process. Though the efficiency 
for selecting an inclusive v^ — GC sample and the 



resulting sample purity are not strongly depen- 
dent on the re-scattering model, we find that 
the purity of a quasi-elastic sub-sample varies by 
~ 10% when we compare simulations with and 
without rescattering. We expect that efforts to 
select other exclusive final states (such as coher- 
ent 7r production) or to discriminate v e from NC, 
will be similarly afflicted. 

The model employed here is under active devel- 
opment intended to widen its scope and increase 
its sophistication. The addition of baryon re- 
scattering seems quite important as does a more 
detailed treatment of the energy lost in inelastic 
collisions. Though the model is anchored to v + d 
and v + Ne data, further benchmarking, either 
against additional data or theoretical predictions, 
would be very valuable. For example, the author 
would be interested in independent predictions of 
Figures [21 and 
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Figure 2. The probability, according to NEUGEN, 
that a pion produced in a neutrino interaction on 
56 Fe undergoes at least one inelastic interaction 
(full circles) or is absorbed (open circles) before 
exiting the nucleus. 
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Figure 3. The average momentum lost, due to 
intranuclear re-scattering, by pions produced in 
neutrino interactions on 56 Fe. The curves are the 
result of NEUGEN simulations. The full circles cor- 
respond to all pions exiting the nucleus while the 
open circles correspond to only those pions that 
scattered while exiting. 
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Figure 4. The variation in the calorimeter's re- 
sponse as a function of the energy transfer v 
and under different final state interaction sce- 
narios. Zero on the vertical scale corresponds 
to the 7T — > npnp model of pion absorption and 
to=0.52 fm/c. 



Figure 5. The reconstructed energy spectrum of 
events selected as v^ — CC under three different 
final state interactions scenarios. The events were 
generated with NEUGEN and simulated in a full 
GEANT simulation of the MINOS Near Detector. 
The underlying spectrum is that expected from 
the NuMI low-energy beam configuration. 
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